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ABSTRACT 


hysical  acoustic  sources  that  radiate  transient  waveforms  represent  a 
large  class  of  noise  generators.  Active  transmissions,  flow-noise,  and 
helicopter-radiated  noise  under  certain  flight  situations  depending  on  the 
source-observer  orientation  are  members  of  this  class.  The  objective  of 
this  paper  is  to  present  a  signal  processing  methodology  to  extract  time 
domain  recurring  waveforms  from  data.  The  recurring  waveforms,  on  the  one 
hand,  may  represent  a  source  of  interference  that  is  to  be  eliminated  as 
much  as  possible.  But,  on  the  other  hand,  these  recurring  waveforms  may 
represent  a  desired  signal  that  is  to  be  extracted  from  a  background  of 
undesired  noise.  Both  cases  are  shown  in  the  paper  to  be  different  aspects 
of  the  same  problem.  Initially,  the  methodology  is  based  on  a  physical 
property  of  equally  spaced  and  identical  recurring  waveforms.  These  results 
are  then  generalized,  because  in  practice,  either  by  design  or  by  other 
causes,  the  recurring  waveforms  may  not  be  equally  spaced  and  identical. 

The  method  is  then  applied  to  extract  equally  spaced  and  nearly  identical 
waveforms  which  were  produced  by  an  approaching  helicopter. 

\ 

ADMINISTRATIVE  INFORMATION 

This  technical  memorandum  was  prepared  under  NUSC  Project  No.  A75031, 
"Multidimensional  Signal  Processing,"  Principal  Investigator  Or.  Roger  F. 
Dwyer  (Code  3314),  Program  Manager  Or.  N.  Gerr,  ONR  (Code  411). 


The  author  of  this  document  is  located  at  the 
New  London  Laboratory,  Naval  Underwater  Systems  Center,A 

New  London,  Connecticut  06320.  sM-'i  , 

'  LVi  ■ 

I 

i  '.i' 4  ' . 


i  J_'  ■ 


V 

3/4 

Reverse  Blank 


TM  No.  861043 

INTRODUCTION 

The  data  to  be  utilized,  later  In  the  paper,  as  an  application  to  the 
signal  processing  methodology  Is  from  a  helicopter.  Therefore,  a  short 
discussion  of  aerodynamic  noise  generation  appropriate  for  the  application 
Is  presented  Initially.  A  more  general  discussion,  then  the  one  presented, 
can  be  obtained  from  the  references  cited.  The  discussion  presented, 
however,  will  lead  to  the  signal  processing  methodology  which  Is  the  main 
topic  of  the  paper.  In  addition,  the  signal  processing  methodology  can  be 
applied  to  active  transmissions  and  flow-noise  problems. 

The  theory  of  aerodynamic  noise  generation  Is  based  on  Llghthill's  [1] 
formulation  of  the  wave  equation 
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> 


These  sources  on  the  right  of  equation  1  are  described  In  more  detail 
by  Lowson  and  applied  to  helicopter  noise  [2].  Lowson  points  out  that  the 
first  term  represents  a  simple  source  produced  by  mass  Introduction  at  a 
point  which  produces  a  monopole  source.  The  second  term  represents  sound 
produced  by  forces  acting  on  the  medium  and  can  be  described  as  a  dipole 
source.  Whereas,  the  last  term  Is  due  to  acoustic  stresses  and  they  radiate 
as  a  quadrupole  source. 


For  a  rigid  surface  In  motion  Ffowcs-Wllllams  and  Hawkins  [3]  obtained 
the  exact  solution  of  equation  (1)  as  follows: 


(2) 


where  Is  the  outward  normal  force  per  unit  area  exerted  on  the 

fluid  by  the  surface  of  the  body  S,  u^  Is  the  outward  normal  component  of 
the  surface  convection  velocity  U,  and  Is  the  uniform  density  of  the 
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This  equation  was  evaluated  for  rotors  operating  at  supersonic  tip 
speeds  by  Hawkings  and  Lowson  [4].  In  this  case  they  assumed  the  blades 
were  thin  and  therefore  the  last  term  of  equation  (2)  could  be  neglected. 

Since  the  acoustic  radiated  noise  from  the  rotor  is  periodic  they 
decomposed  it  into  Its  harmonic  components  using  Fourier  analysis.  Although 
more  general,  their  results  also  agreed  with  results  of  Gutin  [4]. 

The  sources  of  aerodynamic  radiated  noise  from  a  helicopter  can  be 
grouped  into  three  categories  [5,  6,  7]:  1)  Rotational  noise,  or  Gutin 

noise,  which  produces  discrete  frequencies  at  multiples  of  the  blade  passage 
frequency,  2)  Broadband  noise  due  to  blade-turbulance  interaction  and  3) 
impulsive  noise  due  to  blade-vortex  interactions  and  high  tip  speeds  which 
also  produce  discrete  frequencies,  at  multiples  of  the  blade  passage 
frequency,  but  at  significantly  higher  levels  than  for  Gutin  noise. 

In  equation  (2)  the  square  brackets  mean  that  the  contents  are  to  be 
evaluated  at  the  retarded  time,  T=t  -r/c,  where  r  =|X  -  Y|1s  the  distance 
from  the  source  Y  to  the  observer.  The  other  term,  1/|1-M^|,  is  the 
doppler  factor  where  is  the  component  Mach  number  in  the  direction  of 
the  observer. 

An  important  physical  insight  to  understanding  the  mechanism  of 
helicopter  impulsive  noise  was  contributed  by  Schmitz  and  Yu  [8].  They 
showed  that  the  doppler  factor  gave  a  large  pressure  amplification  as  a 
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function  of  the  azimuthal  source  position  of  the  helicopter  blade.  Their 
description  accounts  for  both  the  recurring  impulses  and  directionality 
associated  with  helicopter-radiated  impulsive  noise.  Essentially  they 
showed  that  for  an  approaching  helicopter  impulses  are  heard  in  and  near  the 
plane  of  the  main  rotor  and  disappear  entirely  directly  beneath  the  rotor. 
They  attributed  the  mechanism  of  impulsive  noise  to  sound  arriving  at  the 
observer's  position  in-phase  due  to  the  sound  emitted  when  the  approaching 
blade's  azimuthal  position  was  almost  perpendicular  to  the  direction  of  the 
observer.  For  large  component  Mach  numbers  (approaching  one)  in  the 
observer's  direction  Impulsive  noise  would  be  heard  and  would  recur  at  the 
blade  passage  frequency. 

The  time-domain  theoretical  results  of  reference  [8]  were  compared  with 
in-flight  tests  of  an  UH-1H  helicopter  [9]  and  with  wind  tunnel  tests  [10]. 
The  theoretical  results  were  found  to  be  in  good  agreement  with  the  real 
data  except  that  the  large  negative  pressure  peaks  associated  with  Impulsive 
noise  were  always  under  predicted. 

The  objective  of  this  paper  is  to  present  a  signal  processing  method  to 
extract  recurring  waveforms  from  a  system  so  as  to  enhance  a  noisy 
measurement  or  to  eliminate  them  as  a  source  of  interference. 

Field  measurements  may  be  contaminated  by  other  noise  sources  or  the 
signal  may  be  weak  if  measured  at  long  ranges.  In  these  cases  it  would  be 
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advantageous  to  be  able  to  enhance  the  time-domain  data  and  Improve 
interpretation  and  comparison  with  other  data. 


RECURRING  TRANSIENTS 


The  physical  principle  of  extracting  recurring  waveforms  from  data  is 
contingent  on  obtaining  discrete  frequencies  from  its  spectrum.  Therefore, 
the  methodology  is  based  on  frequency-domain  processing  [11].  This  method 
relies  on  obtaining  many  equally  spaced  and  identical  waveforms  over  an 
Interval  of  time.  The  spectrum  of  this  data  over  the  interval  approaches  a 
harmonically  related  line  spectrum  as  the  number  of  waveforms  increase.  The 
discrete  frequencies  contain  all  the  useful  information  of  the  recurring 
waveforms  over  the  interval.  An  enhancement  occurs  if  the  data  of  interest 
are  composed  of  a  relatively  few  line  components  distributed  over  a  larger 
band.  Therefore,  this  method  enhances  the  time-domain  data.  In  addition, 
in  practice,  this  procedure  will  always  be  an  approximation.  Nevertheless, 
an  example  of  a  helicopter  producing  recurring  waveforms  presented  later 
will  show  that  the  method  is  useful  even  for  a  short  time  interval  of  data. 
But,  once  the  enhanced  time  domain  data  is  available,  subsequent  frequency- 
domain  processing  may  be  applied. 

Consider  an  acoustic  waveform  recurring  every  t'  seconds.  The  pressure 
received  at  a  sensor  at  time  t  can  be  expressed  by  the  following  equation 
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go 

Pr(t)  P(t-nt' )  . 

n=-oo 


or  in  the  frequency  domain 


Pp(w) 


4(w-n(i)^) 


n=-oo 


(3) 


where  =  2*/t',  and  P(w)  is  the  Fourier  transform  of  the  transient 
pressure  p(t) . 

Equation  (3)  implies  that  for  a  recurring  waveform,  p(t),  its 
associated  frequencies  are  all  discrete  and  are  harmonically  related  to 

For  the  case  of  helicopter-radiated  noise  is  the  blade  passage 
frequency  [6].  But  in  practice  we  may  not  have  an  infinite  number  of 
recurring  waveforms.  However,  for  a  stationary  source  [6,  7,  12]  this  can 
be  approximated  to  a  high  degree.  In  the  experiments  of  Schmitz  and  Boxwell 
[9]  they  were  able  to  obtain  a  large  number  of  recurring  waveformsfor  a 
source  in  motion.  In  this  case,  the  frequencies  were  doppler  shifted. 


For  the  data  to  be  presented  later  the  source  was  approaching  a  fixed 
observer  and  therefore  the  radiated  wavefomfewere  nonstationary  over  time 
with  respect  to  the  fixed  observer.  Therefore,  only  over  a  short  interval 
of  time  could  the  data  be  assumed  stationary. 

In  the  rest  of  the  paper  the  data  are  assumed  to  be  filtered  to  the 
desired  passband  and  appropriately  A/0  converted. 

Let  the  received  discrete  data  be  given  by  the  expression 

x(i)  h[i-(igm')]  (4) 

k=0 

where  i-0,1 ,2, . . . ,N-1  are  the  discrete  temporal  samples  and  w  is  the  number 
of  waveforms  occurring  over  the  interval  of  N  uniformly  spaced  samples.  For 
convenience  the  sampling  interval  is  assumed  to  be  equal  to  one.  The  other 
parameters  are  as  follows; 

ig  -  The  sample  on  which  the  first  waveform  starts. 

t'  -  Temporal  separation  of  waveforms,  given  as  an  integral 

number  of  samples. 

h  -  Basic  shape  of  waveform. 

A.  '  Amplitude  of  k-th  waveform. 
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Since  the  waveforms  are  constrained  to  occur  over  a  finite  Interval  of 
N  samples  the  following  Inequality  must  hold 

0  <  kf  <  N-l-1g-a 


where  A  Is  the  width  of  the  waveform. 


Transforming  x(i)  to  the  frequency  domain  we  obtain  the  spectrum 


N-1 

X(K)  =  (1/jN)  ^  x(1)  ex(i)  (5) 

1=0 


where  lf=0,l,2 . N-1  and  we  have  defined  the  expression  ex(  )  =  e 


The  spectrum  of  equation  (4)  Is  given  by 

r 


w-l 

X(K)  =  (WyjN)  ^  W(K,k)  ex(i  +kt') 
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where 


+kt') 

A.h(r)  ex(r)  (6) 

+kt’) 


However,  for  the  class  of  functions,  h{r),  being  considered  here 


N-l-(i 

W(K,k)  =  3 


W(K,k)  =  A  W(K). 


For  example,  let 


h(r)  =  -  (1/4)  -  (1/2)  6(r)  -  (1/4)  Mr-t^) 

+  (1/2)  4(r+t2-t^)  4(r-t^)  +  (1/2)  4(r-t^-t^) 

-  (1/4)  4(r^-t2-2t^)  -  (1/2)  4(r-2t^)  -  (1/4)  4(r-t2-2t^) 

where  <$(  )  is  the  Kronecker  delta  function  and  the  following  inequalities 
hold: 


0  <  <  21^  <  t,  <  Jt,  <  a  . 

0  <  a  <■  i  4-  kt'  <  N-1 

—  n  — 


,  k=0,l ,2, . . , ,w-l . 


TM  No.  861043 


The  function,  h(r).  Is  plotted  In  figure  1  for  three  recurring  impulse 
functions.  We  also  constrained  h(r)  to  have  an  average  value  equal  to  zero 
over  the  interval.  This  will  insure  that  its  power  spectrum  is  zero  at  K=0. 

From  equation  (6)  we  obtain  the  expression. 


W(K)  =  [Ucos(2irKt2/N)]  [1 -COS(  2irKt^ /N)  ]  ex(t^) 


and  the  spectrum  for  this  example  is  given  by 


X(K)  =  (Aw/yN)  [l+COS(2irKt2/N)]  [1 -COS(2irKt^N)  ]  ex(t.|) 


s1n(irwKt'/N) 
wsin(irKt '  /N) 


(7) 


The  peaks  of  equation  (7)  can  be  conveniently  found  from  its, 
unaveraged,  magnitude  squared  or,  power  spectrum,  which  is  given  by 

|X(K)|^  =  (aV/N)  [l4-cos(2irKt2/Nl^  [1 -COS(2irKtyN) 

sin(irwlct'/N) 

■  w  sin(irKt'/N) 
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where  the  expression  [1+cos(2irKt2/N)]^  [l-cos(2irKt.|/N]^ 
controls  the  overall  shape  by  the  locations  of  Its  zeros.  Whereas, 


sinU  wKf/Nl 
w  sin(irKt'/N) 


give  the  locations  of  the  peaks  of  the  power  spectrum  from  the  solutions  to 
the  equation 


s1n(irKt'/N)  *  0 


which  are  given  by 


K  =  (N/f) 


for,  0  <  k  <  N-1  and  I  =  0,  +1,  +2,  ...  . 

—  —  0  ~  “ 
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Therefore,  t'  controls  both  the  number  and  locations  of  the  discrete 
components.  As  w  equation  (9)  will  approach  a  unit  delta  function  so 
that  the  spectrum  becomes  a  purely  line  spectrum. 

The  power  spectrum  of  equation  (8)  is  plotted  in  figure  2  for  set  of 
realistic  parameters  which  correspond  with  the  helicopter  data  that  will  be 
presented  later.  The  parameter  set  is  as  follows;  N  =  1024,  w  =  12,  A  =  .3 
Volts,  T^  =  1  msec,  =  A  msec,  and  t'  =  83  msec. 

EXTRACTION 

The  previous  example  showed  that  a  recurring  process  produced  a  power 
spectrum  with  sharp  peaks.  Therefore,  most  of  the  information  was  contained 
at  these  peaks.  This  phenomenon  is  the  basis  for  the  extraction  method. 

Rewriting  equation  (7)  in  the  following  way 


reveals  the  extraction  procedure.  The  term 


sinlirwICt'/N) 

wsinlirlCt'/N) 


16 


TM  No.  861043 


represents  the  locations  of  the  peaks.  These  solutions  form  a  solution 
space  denoted  by  S.  The  discrete  frequencies  representing  the  solution 
space  will  be  denoted  by  the  expression  KcS. 

The  extracted  data  are  found  by  transforming  X(K),  for  all  KcS,  into 
the  time-domain  by  applying  an  inverse  discrete  Fourier  transform  (IDFT)  as 
follows: 


^  X(K)ex(-i) 

KcS 

and  from  equation  (10), 

x^(i)  =  (Aw/N)  ^  W(lf)  ex|-[i-  (ig+’i^  t')]f  (11) 

KeS 

In  practice  the  inverse  fast  Fourier  transform  is  utilized  so  that  the 


method  is  highly  efficient. 
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Now  the  resultant  data,  left  after  the  extraction  process.  Is  obtained 
from  the  equation 


Xp  (1)  =  x(1)  -  x^(i) 


(12) 


where. 


x^(1)  =  (Aw/N) 


E 


W(l^)  exl-[i-(1  + 


w-1 

2 


Ul^t'  /Nil 


(irlft'/N)J 


SO  that  the  original  data  Is  the  sum  of  the  extracted  and  resultant 
components, 


x(i)  =  x^(i)  +  Xjj(i). 


These  Ideas  can  be  discussed  in  terms  of  matrices  and  vectors.  This  is 
called  the  matrix  method. 
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Let 


X  =  (x(0),  x(l) . x(N-l))^ 


be  a  vector  of  time  domain  data  samples.  Here  T  represents  transpose. 


Similarly,  the  frequency  domain  data  is  given  by  the  vector 


X  =  (X(0),  X{1),  ....  X(N-1))^ 


where  the  two  vectors  are  related  through  the  transformation 


X  =  F  X 


and  F  is  a  (NxN)  transform  matrix  representing  the  DFT  operation  as  follows: 
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F 


1  1 
1  exp(-j2ir/N) 

1  exp(-j2ir2/N) 


1 

exp{-j2ir2/N) 

exp(-j2»4/N) 


1  exp(-j2ir(N-2)/N)  exp( -j2ir2(N-2) /N  .  .  . 
1  exp(-j2ir(N-l)/N)  exp(-j2ir2(N-1)/N)  .  .  . 


and  Its  Inverse  exists  and  Is  denoted  by 


=  (1/N) 


1/2 


1  1  1 
1  exp(j2ir/N)  exp(j2ir2/N) 
1  exp(j2ir2/N)  exp(j2^4/N) 


1  exp(j2tr(N-2)/N)  exp( j2«2(N-2)/N)  .  . 
_1  exp(j2ir(N-l)/N)  exp(j2#2(N-l)/N)  .  . 


1 

exp(-j2ir(N-l)/N) 

exp(-j2ir2(N-l)/N) 


exp(-j2ir(N-2)(N-l)/N) 

exp(-j2ir(N-l)(N-l)/N) 


1 

.  exp(j2ir(N-l)/N) 

.  exp{j2ir2{N-l)/N) 


.  exp(j2ir(N-2)(N-l)/N) 

.  exp(j2ir(N-l)(N-l)/N) 
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We  will  now  define  a  diagonal  (NxN)  operator  matrix  Q  which  is  applied 
in  the  frequency  domain. 


The  frequency  domain  resultant  data  is  then  represented  in  the  form 


^  -  [F-QF]  X 


Applying  the  inverse  transformation  we  obtain 


x_  =  X  -  F  QF  X. 


If  0=1  (identify  matrix),  then,  =  0  (zero  vector). 
Whereas,  if  Q  =  0  (zero  matrix),  then,  Xj^  =  x. 

Therefore  we  can  choose  0  to  obtain  the  desired  result. 


In  addition. 
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%  =  F'^QP  X 


is  the  extracted  data  vector 

In  applications  the  elements  of  diagonal  matrix  Q  are  chosen  to 
represent  the  locations  where  the  discrete  components  occur  in  the  frequency 
domain  i.e.,  KcS. 

For  example,  let  N=1024,  and  let  there  be  four  exact  solutions  located 
at  K=0.  256,  512,  and  768, 

Then  Q,  =  diag  1.j=0,l .  1023] 

where. 


1,  at  1=j  =  0,  256,  512,  and  768 
0,  otherwise 


Therefore, 
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*i»here, 

rioooioooi  .  . 

I 

0  1  0  0  0  1  0  0  0  1  . 
00100010001 
00010001000 
10001000100 

E  =  . 

0 

0 

jo 

jlOOOlOOOl  ,  . 
0  1  0  0  0  1  0  0  0  1  . 
00100010001 
00010001000 


.0  0  0  1  0  0  0 
.  .0  0  0  1  0  0 

.  .  .0  0  0  1  0 

1  ...  0  0  0  1 
0  1  .  .  .000 


.0  0  0  1  0  0  0 

.  .0  0  0  1  0  0 

.  .  .0  0  0  1  0 

1  .  .  .  0  0  0  1_ 

(1024  X  1024) 


So  the  extracted  data  vector  for  this  example  is  found  from  the 
expression 


Xr  =  (4/N)  E  X  =  X 
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The  extracted  data  vector  was  equal  to  the  original  data  vector  because 
In  this  example  all  the  information  was  located  at  discrete  frequencies. 

Now  if  we  add  noise  of  the  form  n  then 


2  =  X  +  n 


and 


2^  =  (4/N)  E  X  +  (4/N)  E  n 

=  X  f  (4/N)  E  n 


(13) 


which  reduces  to,  2^  ~  x.  Since  the  noise  is  random,  the  noise 
vector  will  not  match  the  transformation  E  and  consequently  the  second  term 
on  the  right  side  of  equation  (13)  will  be  small. 

The  resultant  data  vector  is  obtained  from  the  expression. 
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Adding  noise  again  we  obtain 


Zp  =  n  -  (4/N)  E  n  ~  n 


(14) 


Therefore,  for  this  example,  the  extraction  method  is  nearly  perfect 
even  with  additive  noise.  In  practice  the  matrix  method  will  be  limited  by 
the  ability  to  choose  the  appropriate  operator  matrix  Q.  In  this  regard  we 
have  shown  the  physical  connection  between  the  peaks  in  the  spectrum  (or 
solutions  of  equation  (9))  and  the  operator  matrix.  But  the  matrix  method 
is  a  purely  mathematical  procedure.  Therefore,  we  can  generalize  this 
procedure  to  handle  other  situations  that  can  occur  in  the  time  domain  but 
do  not  have  an  obvious  physical  connection  in  the  frequency  domain. 

Therefore,  let  E^  be  a  (NxN)  matrix  with  elements 


[<3^y  i.  j=0,l,2 . N-1]. 


Then, 


x_  =  E  X  and  x„  =  x  -  E  x 
-E  g  -  -R  -  g  - 
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are  the  extracted  and  resultant  vectors  for  the  generalized  transformation 

matrix  E  . 

9 

In  this  generalized  formulation  x  need  not  represent  a  vector  of 
uniformly  spaced  and  identical  recurring  waveforms.  Rather,  the  recurring 
waveforms  can  have  unequal  separations  with  arbitrary  amplitudes.  In  fact, 
the  shapes  need  not  be  the  same  over  the  interval.  However,  the  elements  of 
the  generalized  transformation  matrix  must  contain  this  knowledge  in  order 
to  match  the  data  vector.  This  information  may  be  obtained,  perhaps,  from  a 
learning  procedure  or,  on  the  other  hand,  it  may  be  known  a  priori  as  in  a 
coded  message. 

In  the  next  section  we  shall  consider  a  segment  of  real  data  that 
contain  recurring  waveforms  produced  in  a  system  by  an  approaching 
helicopter. 


RESULTS 


The  details  of  the  data  acquisition  were  discussed  in  reference  11. 
Essentially,  data  from  an  approaching  UH-1  helicopter  were  recorded  from  a 
microphone  mounted  2  meters  above  the  ground.  A  one  second  segment  of 
far-field  data  were  analyzed  after  being  bandpassed  filtered  and  sampled  at 
1024  samples  per  second  using  a  12  bit  A/0  converter.  The  band  was  from 
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The,  unaveraged,  power  spectrum  estimate  of  this  one  second  interval  of 
data  is  shown  in  figure  3.  A  1024-point  FFT  was  used  to  obtain  the 
frequency  domain  data.  Therefore,  the  frequency  resolution  was  1  Hz.  The 
power  spectrum  in  figure  3  is  characteristic  of  helicopter-radiated  noise  as 
discussed  in  reference  5.  Here,  however,  only  the  main  rotor  frequency 
components  are  discernable  in  the  figure.  The  peaks  are  indicated  by 
arrows.  Apparently,  this  was  due  to  the  range  of  the  helicopter. 

Subsequent,  closer-in,  segments  revealed  contributions  from  both  main  and 
tail  rotors. 

These  eight  frequency  components  indicated  in  the  figure  by  arrows  will 
be  utilized  to  obtain  the  extracted  and  resultant  data  components. 

Actually,  however,  the  adjacent  samples  will  be  included.  This  gives  a 
total  of  48,  out  of  the  possible  1024  samples,  incorporated  in  the  extracted 
components. 

During  the  approach  of  the  helicopter  impulsive  noise  could  be  clearly 
heard.  These  recurring  waveforms  are  shown  in  the  top  graph  of  figure  4. 
Over  the  one  second  interval  there  are  12  recurring  waveforms.  Whereas,  the 
time  interval  between  waveforms  is  82  msec.  This  corresponds  to, 
approximately,  a  blade  passage  frequency  of  12  Hz.  The  indicated 
frequencies  of  figure  3  are  all  multiples  of  12  Hz.  But  due  to  the  bandpass 
filter  employed,  the  multiples  below  50  Hz  are  not  discernable  and, 
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therefore,  not  included  in  the  extraction  procedure.  This  means  that  the 
fundamental  frequency  is  not  needed  to  reconstruct  the  proper  time 
occurrences  of  the  waveforms.  The  reason  is  that  this  information  is 
already  contained  in  the  eight  discernable  frequency  multiples.  However,  to 
reconstruct  the  individual  waveform  exactly  all  the  frequencies  associated 
with  the  pulse  must  be  extracted. 

The  middle  graph  of  figure  4  represents  the  resultant  data  after  the 
eight  indicated  frequencies  are  extracted.  Whereas,  the  extracted  data  are 
shown  in  the  bottom  graph.  From  these  plots  it  is  clear  that  the  eight 
indicated  frequencies  represent  the  original  waveforms  to  an  extent,  tiut 
the  result  are  only  an  approximation.  The  time  occurrences  are  maintained 
in  the  extracted  data  but  the  shape  is  distorted.  Nevertheless,  depending 
on  the  application,  these  results  could  be  beneficial.  For  example,  in  a 
noisy  environment  the  time  domain  data  may  be  obscured.  As  long  as  the 
frequency  components  are  still  discernable  the  extracted  data  would  give  an 
enhancement.  As  was  shown  in  the  theoretical  discussing,  extending  the 
interval  and  including  more  waveforms  would  increase  the  discernability  of 
the  frequency  components. 

The  possible  enhancement  is  demonstrated  in  figure  5.  The  original 
data  with  additive  computer  noise,  at  a  -7  dB  signal -to-noise  level 
(measured  in  the  time  domain),  is  shown  is  the  top  graph.  With  this 
signal-to-noise  level  the  eight  frequency  components  were  still 
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discernable.  As  in  figure  4,  the  middle  graph  is  the  resultant  data  after 
the  eight  frequency  components  are  extracted.  The  bottom  graph  representing 
the  extracted  data  clearly  shows  the  time  domain  enhancement. 

The  sensitivity  of  the  extraction  procedure  to  random  Gaussian  noise  is 
demonstrated  in  figure  6.  Here  an  analog  Gaussian  source  was  processed  in  a 
similar  manner  to  the  original  helicopter  data.  As  in  the  previous  figures, 
the  top  graph  is  the  original  Gaussian  noise  data.  The  middle  and  bottom 
graphs  are,  respectively,  the  resultant  and  extracted  data.  These  results 
also  correspond  with  the  theoretical  predictions  of  sensitivity  to  random 
noise  given  by  the  matrix  method. 


SUMMARY  AND  CONCLUSIONS 


In  this  paper  a  signal  processing  methodology  was  presented  to  extract 
recurring  waveforms  from  data.  The  procedure  was  based  on  the  principle 
that  information  associated  with  recurring  waveforms  is  contained  in 
discrete  frequencies.  Therefore,  it  was  assumed,  initially,  that  the 
recurring  waveforms  were  uniformly  spaced  and  identical.  Under  this 
assumption  a  physical  connection  between  time  and  frequency  domains  was 
demonstrated.  However,  from  the  time  domain  development  a  purely 
mathematical  generalization  of  the  procedure  was  revealed  based  on  a 
transformation  matrix.  This  allowed  a  relaxation  of  the  constraints 
originally  required  for  the  extraction  procedure. 
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The  signal  processing  procedure  was  demonstrated  with  helicopter- 
radiated  noise  data.  It  was  shown  that  the  extraction  procedure  could  be 
useful  to  either  remove  recurring  waveforms  from  data  of  interest  or,  on  the 
other  hand,  to  enhance  a  noisy  time  domain  signal.  These  results  are 
currently  being  generalized  for  an  array  of  sensors. 
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